#ifndef _MPoint
#define _MPoint
//
//-
// ==========================================================================
// Copyright (C) 1995 - 2006 Autodesk, Inc., and/or its licensors.  All 
// rights reserved.
// 
// The coded instructions, statements, computer programs, and/or related 
// material (collectively the "Data") in these files contain unpublished 
// information proprietary to Autodesk, Inc. ("Autodesk") and/or its 
// licensors,  which is protected by U.S. and Canadian federal copyright law 
// and by international treaties.
// 
// The Data may not be disclosed or distributed to third parties or be 
// copied or duplicated, in whole or in part, without the prior written 
// consent of Autodesk.
// 
// The copyright notices in the Software and this entire statement, 
// including the above license grant, this restriction and the following 
// disclaimer, must be included in all copies of the Software, in whole 
// or in part, and all derivative works of the Software, unless such copies 
// or derivative works are solely in the form of machine-executable object 
// code generated by a source language processor.

// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND. 
// AUTODESK DOES NOT MAKE AND HEREBY DISCLAIMS ANY EXPRESS OR IMPLIED 
// WARRANTIES INCLUDING, BUT NOT LIMITED TO, THE WARRANTIES OF 
// NON-INFRINGEMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE, 
// OR ARISING FROM A COURSE OF DEALING, USAGE, OR TRADE PRACTICE. IN NO 
// EVENT WILL AUTODESK AND/OR ITS LICENSORS BE LIABLE FOR ANY LOST 
// REVENUES, DATA, OR PROFITS, OR SPECIAL, DIRECT, INDIRECT, OR 
// CONSEQUENTIAL DAMAGES, EVEN IF AUTODESK AND/OR ITS LICENSORS HAS 
// BEEN ADVISED OF THE POSSIBILITY OR PROBABILITY OF SUCH DAMAGES. 
// ==========================================================================
//+
//
// CLASS:    MPoint
//
// *****************************************************************************
//
// CLASS DESCRIPTION (MPoint)
//
//	This class provides an implementation of a point.  Numerous convienence
//	operators are provided to help with the manipulation of points. This
//	includes operators that work with the MVector and MMatrix classes.
//
// *****************************************************************************

#if defined __cplusplus

// *****************************************************************************

// INCLUDED HEADER FILES



#include <maya/MTypes.h>
#include <maya/MStatus.h>
#include <maya/MVector.h>

// *****************************************************************************

// DECLARATIONS

class MMatrix;
class MVector;
class MFloatPoint;
class MFloatVector;

#define MPoint_kTol	1.0e-10

// *****************************************************************************

// CLASS DECLARATION (MPoint)

/// Implementation of a point. (OpenMaya) (OpenMaya.py)
/**
  This class implements the Maya representation of a point
*/
#ifdef _WIN32
#pragma warning(disable: 4522)
#endif // _WIN32

class OPENMAYA_EXPORT MPoint  
{
public:
	///
					MPoint();	// defaults to 0,0,0,1
	///
					MPoint( const MPoint & srcpt );  
	///
					MPoint( const MFloatPoint & srcpt );  
	///
					MPoint( const MVector & src );  
	///
					MPoint( const MFloatVector & src );  
	///
					MPoint( double xx, double yy,
							double zz = 0.0, double ww = 1.0 );
	///
					MPoint( const double[4] );
	///
					MPoint( const float[4] );
	///
					~MPoint();
	///
	MStatus			get( double[4] ) const;
	///
	MStatus			get( float[4] ) const;
	///
	double  		operator[](unsigned int i) const;
	///
	double  		operator()(unsigned int i) const;
	///
	MPoint & 		operator=( const MPoint & src );
	///
	MVector   		operator-( const MPoint & other ) const;
	///
	MPoint   		operator+( const MVector & other ) const;
	///
	MPoint   		operator-( const MVector & other ) const;
	///
	MPoint & 		operator+=( const MVector & vector );
	///
	MPoint & 		operator-=( const MVector & vector );
	///
	MPoint			operator*(const double scale) const;
	///
	MPoint			operator/(const double scale) const;
	///
	MPoint   		operator*(const MMatrix &) const;
	///
	MPoint & 		operator*=(const MMatrix &);
	///
	bool           	operator==( const MPoint & other ) const;
	///
	bool           	operator!=( const MPoint & other ) const;
	///
	MPoint & 		cartesianize();
	///
	MPoint & 		rationalize();
	///
	MPoint & 		homogenize();
	///
	double     		distanceTo( const MPoint & other ) const;
	///
	bool           	isEquivalent( const MPoint & other,
								  double tolerance = MPoint_kTol) const;

BEGIN_NO_SCRIPT_SUPPORT:

	///	NO SCRIPT SUPPORT
	double &        operator[](unsigned int i);

	///	NO SCRIPT SUPPORT
	double &        operator()(unsigned int i);

	///	NO SCRIPT SUPPORT
	friend OPENMAYA_EXPORT MPoint operator*(const MMatrix &, const MPoint &);
	///	NO SCRIPT SUPPORT
	friend OPENMAYA_EXPORT std::ostream& operator<<(std::ostream& os, const MPoint& p);

END_NO_SCRIPT_SUPPORT:

	///
	static const MPoint origin;
public:
	/// the x component of the point
	double			x;
	/// the y component of the point
	double			y;
	/// the z component of the point
	double			z;
	/// the w component of the point
	double			w;

protected:
// No protected members

private:

	static const char* className();
};

#ifdef _WIN32
#pragma warning(default: 4522)
#endif // _WIN32

inline MPoint::MPoint()
 	: x(0.0)
	, y(0.0)
	, z(0.0)
	, w(1.0)
{
}

inline MPoint::MPoint(const MPoint& srcpt)
 	: x(srcpt.x)
	, y(srcpt.y)
	, z(srcpt.z)
	, w(srcpt.w)
{
}

inline MPoint::MPoint(const double d[4] )
 	: x(d[0])
	, y(d[1])
	, z(d[2])
	, w(d[3])
{
}

inline MPoint::MPoint(const float d[4] )
 	: x(d[0])
	, y(d[1])
	, z(d[2])
	, w(d[3])
{
}

inline MPoint::MPoint(double xx, double yy, double zz, double ww)
 	: x(xx)
	, y(yy)	
	, z(zz)	
	, w(ww)
{
}

inline MPoint::~MPoint()
{
}

inline MPoint& MPoint::operator=(const MPoint& src)
{
	x = src.x;
	y = src.y;
	z = src.z;
	w = src.w;
	return (*this);
}

inline double& MPoint::operator() (unsigned int i)
{
	switch( i )
	{
		default:
		case 0: return x;
		case 1: return y;
		case 2: return z;
		case 3: return w;
	}
}

inline double MPoint::operator() (unsigned int i) const
{
	switch( i )
	{
		default:
		case 0: return x;
		case 1: return y;
		case 2: return z;
		case 3: return w;
	}
}

inline double& MPoint::operator[]( unsigned int i )
{
	switch( i )
	{
		default:
		case 0: return x;
		case 1: return y;
		case 2: return z;
		case 3: return w;
	}
}

inline double MPoint::operator[]( unsigned int i ) const
{
	switch( i )
	{
		default:
		case 0: return x;
		case 1: return y;
		case 2: return z;
		case 3: return w;
	}
}

inline MStatus MPoint::get( float dest[4] ) const
{
	if( dest != NULL )
	{
		dest[0] = (float) x;
		dest[1] = (float) y;
		dest[2] = (float) z;
		dest[3] = (float) w;
		return MS::kSuccess;
	}
	return MS::kFailure;
}

inline MStatus MPoint::get( double dest[4] ) const
{
	if( dest != NULL )
	{
		dest[0] = x;
		dest[1] = y;
		dest[2] = z;
		dest[3] = w;
		return MS::kSuccess;
	}
	return MS::kFailure;
}

inline MPoint& MPoint::operator+= (const MVector& vector )
{
	x += vector.x;
	y += vector.y;
	z += vector.z;
	return *this;
}

inline MPoint& MPoint::operator-= (const MVector& vector )
{
	x -= vector.x;
	y -= vector.y;
	z -= vector.z;
	return *this;
}

inline MPoint MPoint::operator+ (const MVector& other) const
{
	if(w==1.0) {
		return MPoint(x+other.x, y+other.y, z+other.z);
	} else {
		MPoint p1(*this); p1.cartesianize();
		return MPoint(p1.x+other.x, p1.y+other.y, p1.z+other.z);
	}
}

inline MPoint MPoint::operator- (const MVector& other) const
{
	if(w==1.0) {
		return MPoint(x-other.x, y-other.y, z-other.z);
	} else {
		MPoint p1(*this); p1.cartesianize();
		return MPoint(p1.x-other.x, p1.y-other.y, p1.z-other.z);
	}
}

inline MVector MPoint::operator- (const MPoint& other) const
{
	if(w==1.0 && other.w==1.0) {
		return MVector(x-other.x, y-other.y, z-other.z);
	} else {
		MPoint p1(*this); p1.cartesianize();
		MPoint p2(other); p2.cartesianize();
		return MVector(p1.x-p2.x, p1.y-p2.y, p1.z-p2.z);
	}
}
inline MPoint MPoint::operator* (const double scale) const
{
	return MPoint( x*scale, y*scale, z*scale, w );
}
	
inline MPoint MPoint::operator/ (const double scale) const
{
	return MPoint( x/scale, y/scale, z/scale, w );
}

inline bool MPoint::operator== (const MPoint& other) const
{
	return (x == other.x && y == other.y && z == other.z);
}

inline bool MPoint::operator!= (const MPoint& other) const
{ 
	return !(*this == other); 
}

inline bool MPoint::isEquivalent(const MPoint& other, double tol) const
{
	MPoint diff = *this - other;
	return (diff.x*diff.x + diff.y*diff.y + diff.z*diff.z) < tol*tol;
}

inline double MPoint::distanceTo(const MPoint& other) const
{
	MVector diff = *this - other;
	return diff.length();
}

// *****************************************************************************
#endif /* __cplusplus */
#endif /* _MPoint */
